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We recently started to investigate ray-optical analogs 
of metamaterials. Like metamaterials [1,2], these so- 
called METATOYs ( meta materials for rays) [3] are capa- 
ble of performing positive and negative refraction [4,5]. 
In addition, METATOYs can perform local light-ray 
rotation around the interface normal [6]. 

Here we describe local light-ray rotation around the 
interface normal in terms of Fermat's principle. Fermat's 
principle can be treated as a basic theorem of geomet- 
rical optics, but "is itself only understandable in terms 
of a wave theory" [7]. On the other hand, local light- 
ray rotation around the interface normal has no wave- 
optical analog in the sense that it is not always possible 
to construct a wave in which the phase-front normal - 
the geometrical-optics light-ray direction - has been ro- 
tated as required [3]. What, then, happens if we attempt 
to describe local light-ray rotation around the interface 
normal in terms of Fermat's principle? Here we do just 
that. We find a "natural" formulation of Format's prin- 
ciple in which ray rotation is described by an interface 
between homogeneous media with a complex refractive- 
index ratio. This leads to an extended form of Snell's law 
that uses complex refractive indices. 

It is important to note that the meaning of the imag- 
inary part of the complex refractive index we introduce 
here is different from that commonly used in optics [8] : 
in the former it is associated with ray rotation, in the 
latter with attenuation. 

Format's principle [9] states that a light ray traveling 
between two points takes a path with a stationary optical 
path length, that is, for small variations in the path taken 
the optical path length stays the same. The whole path 
can be calculated using the calculus of variations, but 
simplified calculations can be performed using ordinary 
calculus. 

Figure 1(a) shows a planar interface between two me- 
dia with different refractive indices, n\ and ri2- When 
light travels between two fixed points A and B on cither 
side of the interface, via a point P on the interface, but 



not otherwise fixed, it covers an optical path length 

A = nidi + n 2 d 2 (1) 
= n x <Jx\ + y{ + z\ + n 2 \Jx\ + y% + z\. (2) 

(xi, yi, z{) and {x 2 , y 2 , z 2 ) are the components of the vec- 
tor AP (that is, the vector from A to P) and PB, re- 
spectively, in a Cartesian coordinate system whose (x, y) 
plane coincides with the plane of the interface. We can 




Fig. 1. (Color online) Geometry of refraction at a planar 
interface between two media with different refractive in- 
dices, rii and n 2 . (a) A light ray travels from a point A in 
front of the interface to a point P on the interface, and 
then to a point B behind the interface. The geometrical 
distance between A and P is d\, that between P and B is 
d 2 . (b) The light-ray direction can be represented by two 
angles, 9 and <j>, which respectively represent the angle 
with respect to the interface normal (the z axis) and the 
angle of the projection into the interface plane with re- 
spect to the x axis. Alternatively, the light-ray direction 
can be described by the projection of the normalized di- 
rection vector into the interface plane. With a complex 
plane in the interface plane as shown, this projection can 
then be described by a single complex number, c. 
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use the equations 

X = xi+x 2 , Y = y 1 + y 2 , 



(3) 



which describe the constant separations between the 
fixed points A and B in the x and y directions, to elim- 
inate x 2 and 2/2 from the expression for the path differ- 
ence, Eqn (2). This gives 

A = n 1 yjx* + yl + zf + 

n 2 ^j {X - xtf + (Y - Vl f + zl (4) 

We now find the values X\ and y\ for which the func- 
tion A is stationary. According to Fermat's principle, 
these values then correspond to the point P through 
which the light ray would actually travel. These values 
Xi and y\ have to satisfy the equations 



dA 

dx\ 



0. 



dA 

dyi 



0. 



(5) 



Substitution of the expression for A into these equations 
gives 



n 1 n 2 — = 0, ni n 2 — = 0, 

r\ r 2 ri r 2 



(6) 



1/2 

where rj = + y] + ) (with j = 1, 2). The terms 
can be translated into spherical coordinates </> (the az- 
imuthal angle) and 9 (the angle with the z axis - see Fig. 
1(b)) using the equations 



(7) 



In spherical coordinates, equations (6) are therefore the 
real and imaginary part, respectively, of the equation 

ni sin 9i exp(i0i) = n 2 sin 9 2 e~xp(i(j) 2 ). (8) 

This equation is the basis of the remainder of this pa- 
per. Like in Snell's law, it is not the individual refractive 
indices that matter, but their ratio. As equation (8) is a 
complex equation, it is natural to allow the refractive in- 
dices (and their ratio) to be complex numbers. Equation 
(8) is then an extension of Snell's law: for real refractive- 
index ratios, it describes ordinary refraction according 
to Snell's law, expressing (unlike the Snell's-law formula) 
the fact that the refracted ray lies in the same plane as 
the incident ray and the interface normal; and it leads 
to local light-ray rotation as a natural extension of re- 
fraction with complex refractive-index ratios. We discuss 
these properties below. 

For the following discussion it is useful to visualize the 
extended Snell's law as follows. We interpret the plane 
of the refractive-index interface as a complex plane, cen- 
tered at the point P where the light ray intersects the 
plane, and with the real axis in the x direction and the 
imaginary axis in the y direction, z' is the normal to the 
interface at P. Figure 1(b) shows this coordinate system. 
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Fig. 2. (Color online) Plots of the complex numbers c\ 
and c 2 representing various types of refraction, (a) and 
(b) are examples of standard refraction (n\jn 2 real and 
positive). In (a) n\jn 2 > 1, n\jn 2 < 1 in (b). (c) A 
complex refractive-index ratio n\/n 2 = exp(ia) leads to 
local light-ray rotation through an angle a (here a = 
90°). (d) Light-ray rotation through 180° is equivalent to 
negative refraction with a refractive- index ratio ri\ /n 2 — 
-1. 

We now consider a unit vector in the direction of the 
incident light ray, starting at the origin. We define the 
complex number c as the orthographic projection of this 
unit vector into the complex plane. If we assume that 
light travels in the positive z' direction, then this projec- 
tion uniquely defines the ray direction. For a unit vector 
with spherical-coordinate angles 9 and (f>, c is 



sin 9 exp(i(j 



(9) 



The extended Snell's law, equation (8), can then be 
written in the form 



mci = n 2 c 2 , 



(10) 



where c\ and c 2 is the complex number correspond- 
ing to the incident and refracted ray direction, respec- 
tively. In other words, the product of refractive index 
and the complex number representing the ray direction 
remains constant. We can emphasize the dependence on 
the refractive-index ratio by writing equation (10) in the 
form 

"1 = ^. (11) 
n 2 a 

According to equation (11), the ratio c 2 jc\ is the 
same as the refractive-index ratio n\jn 2 . Therefore, if 
the refractive-index ratio is real then so is the ratio of 
the direction projections. 

This has a simple interpretation. The argument of c\ 
— the spherical-coordinates angle (f>i - defines the plane 
of incidence. Specifically, it describes the angle between 
the plane of incidence and the x axis. This, together with 
the fact that the plane of incidence also contains the z 
axis, completely determines the plane of incidence. A real 
ratio c 2 /c\ means that c 2 lies on the same line through 
the origin as C\ (Fig. 2(a) and (b)), and so the refracted 
ray is also in the plane of incidence. 

Mathematically, it means that (f>i — (f> 2 . This means 
that the exp(i0j) terms (j = 1,2) can be cancelled in 
equation (8), which reduces to Snell's law, 



rii sin 9i = n 2 sin 9 2 . 



(12) 
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Perhaps the simplest examples of complex refractive- 
index ratios arc those of the form 



ni 
n 2 



cxp(ia) 



(13) 



(where a is a real number). According to equation (11), 
this implies that 



— = exp(ia). 

Cl 



(14) 



This means that C2, which characterizes the projection 
of the refracted ray into the interface plane, is rotated 
with respect to c\ through an angle a around the point 
P (Fig. 2(c)). In three dimensions, it means that the 
direction of the refracted ray is that of the incident ray, 
rotated around the z 1 axis through an angle a. This is 
precisely the local light-ray rotation that can be achieved 
with METATOYs [3,6]. 

The case of general complex refractive-index ratios 
ni/n 2 can be approached by writing the left-hand side 
of equation (11) in terms of the modulus and argument 
of this ratio, namely 



(15) 



m _ 


m 












exp 


i arg 


(S)] 


n 2 


n 2 









and the right-hand side in the form 

sin 02 exp(i02) sin 82 exp[i(<^i 
sin 0i exp(i^i) 



sin 9\ exp(i^i) 



sin 82 
sin 9\ 



exp(ia), 
(16) 

which expresses the direction of the refracted light ray's 
projection as that of the incident light ray, rotated 
through an angle a around P. Comparison of the moduli 
of equations (15) and (16) reveals that the change of the 
angle between the ray and the z 1 axis is then given by 
the absolute value of the refractive-index ratio according 
to the equation 



ni 
n 2 



sin 62 
sin 81 ' 



(17) 



which, for real and positive refractive-index ratios, is the 
same as Snell's law (Eqn. (12)). Comparison of the argu- 
ments reveals that the rotation angle a is given by the 
argument of the refractive-index ratio: 



arg 



a. 



(18) 



Now we discuss briefly the case of negative refrac- 
tion [10], for which the refractive-index ratio n\/n 2 is 
real and negative. Negative refraction is fully described 
by Snell's law (and indeed the extended Snell's law), 
from which it then follows that the angle of the ray with 
the z' axis, 8, has to change sign. It can alternatively 
be described by a change of the angle 8 without a sign 
change in combination with a ray rotation around the z' 
axis through 180°. Mathematically, the equivalence be- 
tween negative refraction and ray rotation through 180° 
can be expressed as 



The case n\jri2 = —1 is shown in Fig. 2(d). In fact, a 
Dove-prism-array structure that is ray-optically equiva- 
lent to a refractive-index interface with n\/n 2 = —1 [4] 
is a special case of a Dove-prism-array ray-rotator [6] for 
rotation angle a — 180°. 

Ray rotation is a concept that has no wave-optical 
analog. It is curious that it is possible to describe it - 
in such a natural manner - by using Fermat's principle, 
albeit with a complex refractive-index ratio. More work 
may lead to a deeper understanding. 
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